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Abstract 

A first-principles study of the electronic polarization of BN and A1N nanotubes and their 
graphitic sheets under an external electric field has been performed. We found that the polar- 
ization per atom of zigzag nanotubes increases with decreasing diameter while that of armchair 
nanotubes decreases. The variation of the polarization is related to the exterior angle of the bonds 
around the B or Al atoms rather than that around the N atoms. The increase in the polarization 
of the zigzag nanotubes with decreasing diameter is caused by the large variation of the exterior 
angle when they are wrapped into the tubular form. On the other hand, the decrease in the bond 
length results in the weak polarization of thin armchair nanotubes. 

PACS numbers: 73.22.-f, 77.22.Ej 
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I. INTRODUCTION 

Since the discovery of carbon nanotubes (C NTs)jlj], many experimental and theoreti- 
cal studies have focused on synthesizing nanometer-scale tubular forms of various materials 
and revealing their properties. III-V semiconductors have also attracted considerable atten- 
tion due to their many potential applications in high-power and high-temperature devices, 
because the properties of bulk BN are similar to those of diamond, for example, excel- 
lent thermal conductivity, high chemical resistance, and high melting point. In addition, 
since BN is a piezoelectric material, its tubular structure is a candidate for applications 
in nanometer-scale sensors and actuators. It has so far been reported that single-walled 
BN nanotubes (BN NTs) can be produced by arc discharge [2] and chemical substitution 
react ions [3|]. On the theoretical side, there have been several studies on the polarization of 



C NTs and BN NTs 
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111 ]. For example, Kozinsky and Marzari[l(]] reported 



that the longitudinal polarizability of single-wall C NTs varies as the inverse square of the 
band gap, while the transverse polarizability is proportional to the square of the effective 
radius. Chen et al. [3| proved that a gradual reduction of the band gap of BN NTs occurs 
because of the larger ionicity of the BN bond. In addition, Guo and Lin[s| and Wirtz et a/.j^] 
reported the liner scaling law of the longitudinal polalizabilities with respect to diameter of 
BN NTs. However, there still remains much to be learned about the electronic polarization 
of the NTs composed of III-V elements. 

In this paper, we carry out a first-principles calculation on the l ong itudinal electronic 
polarization of BN and A1N NTs using the Berry phase approach 13|, ll4(. The longitudinal 
electronic polarization per atom of the armchair NTs decreases with decreasing diameter 
while the polarization of zigzag NTs increases and the polarization of the armchair NTs is 
smaller than that of the sheets. This tendency upon polarization cannot be explained by 
the conventional interpretation of the relationship between the energy band gap and the 
polarization [15!], because the band gap of the graphitic sheets is larger than that of the NTs 
and the polarization is expected to increase with the reduction of the band gap. In addition, 
we found that the characteristics of the bond around B or Al atoms have a larger effect on 
the polarization than those of the bond around N atoms. 

The rest of this paper is organized as follows. In Sec. II, we briefly describe the method 
used in this study. Our results are presented and discussed in Sec. III. We summarize our 
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findings in Sec. IV. 



II. COMPUTATIONAL METHODS 



Figure [I] shows the computational models of the III-V NTs and the III-V graphitic sheet. 
The bond length ao is set to 1.43 A for BN tubes and 1.83 A for A1N tubes, which are 
reported by the previous first-principles ca 



real-space finite-difference approach 
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culations[8j, ll6|]. The calculations are performed 
within the local density approximation 17] of density functional theory^, Q] using the 



221 ] and the norm-conserving pseudopotentials 

23, 2J, ]2M. The cutoff 



of Troullier and Martins in the Kleinman-Bylander representation 
energy is set to 110 Ry, which corresponds to a grid spacing h of 0.30 a.u., and a denser grid 
spacing is taken to be h/3 in the vicinity of the nuclei with the augmentation of double-grid 



points 
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221 ] . With this grid spacing, the difference in total energy from that with the grid 



spacing h of 0.15 a.u. is 7 meV/atom. The periodic boundary condition is imposed on the 
supercell in all directions. During the first-principles structural optimization, we relaxed all 
the atoms using a conjugate gradient algorithm, reaching a tolerance in the force of F max < 
lmH/bohr. In the case of the NTs, the length of the supercell along the tube axis is fixed at 
3ao (v^ao) for zigzag (armchair) NTs and a vacuum region of more than 10 A between NTs 
in the neighboring supercells is taken. The integration over the Brillouin zone for the tube 
axis is performed by the equidistant sampling of 10 (18) ^-points for the zigzag (armchair) 
NTs. In the case of the graphitic sheets, the supercell is 3ao X \^3a Q x L Z) where L z is the 
length of the supercell along the z axis, which is taken to be 10 A to avoid the interaction 
between the adjacent sheet. The sampling fc-points over the Brillouin zone is 10 x 18 x 1. 
We have ensured that the enlargement of the supercell and the increase of the number of 
the sampling k points do not change our conclusions significantly. 

An electric field of 0.1 V/ A is applied along the NT axis and along the directions parallel 
to the graphitic sheet. The electronic ground state and the optimized electronic structures 
under external electric fields are determined using the discrete Berry phase scheme proposed 
by Umari and Pasquarello l^j]. The electronic polarization can be expressed as 

P = -2e ^(r,), (1) 

i 

where i is the band index, — e is the electron charge, the factor 2 indicates the spin freedom, 
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and (rj) is the center of the density distribution of the occupied wave functions ipi, an d is 



defined by the following formula 



L 



(rj) = -g-Im (ln(^i|e" 



(2) 



with L being the length of the supercell. The variation of the polarization due to the external 
electric field is defined by AP E = P E — P°, where P E (P°) corresponds to the electronic 
polarization when the external electric field E is present (absent). 



III. RESULTS AND DISCUSSION 

Table [I] lists the electronic polarization per atom of the graphitic sheet. Although the 
graphitic sheets are 120-degrees rotationally symmetrical, the sheets do not exhibit in-plane 
anisotropy. Figure [2] shows the electronic polarization of the NTs. When the atomic ge- 
ometries of the NTs are fully relaxed, the B or Al atoms move toward the central axis of 
the NTs, and the N atoms move in the opposite direction. The bond configuration around 



the B or Al atoms becomes planar as reported by Blase et al. 26j The radial bucklings are 
listed in Table HIl which are agreement with those reported by previous studies 27|, |28( . The 
polarization of the zigzag NTs increases as their diameters become small. In contrast, the 
armchair NTs exhibit the opposite tendency upon polarization; the electronic polarization 
of the armchair NTs becomes smaller than that of the graphitic sheet and decreases as the 
NTs become thin. The different tendency of the polarization with the decrease of the di- 
ameter agrees with the results computed by Guo et a/.jsl, in which the liner scaling law 
of the longitudinal polalizabilities exhibits certain deviations depending on the chiralities. 
Although the different polarization tendency of the BN NTs was found in the previous first- 
principles study [ll!], the origin of the interesting polarization is still unsettled question. It 
is reported that the electronic polarization of conventional III-V compounds increases with 



the reduction of the band gap 151. In addition, all NTs are semiconductors and their band 

n 

gap grows monotonically with increasing diameter [121] . Our results regarding the armchair 
NTs are inconsistent with those predicted from these previous studies because the band gap 
of a III-V graphitic sheet is larger than that of III-V NTs. 

To discuss the different tendency of the polarization between the zigzag and armchair 
NTs in detail, we examine the polarization with respect to the BN bond angle and the cross 
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section of the BN graphitic sheet, which varies during the wrapping of the sheet into the 
tubular form. We first calculate the electronic polarization of the graphitic sheet in which B 
or N atoms are displaced along the direction perpendicular to the graphitic sheet. Schematic 
illustrations of the computational models are shown in Fig. [31 where the neighboring B(N) 
atoms are alternately shifted toward the positive and negative z direction. Figure H] shows 
the polarization as a function of the exterior angle 9 of the BN bonds. The polarization 
increases (decreases) nearly quadratically with the exterior angle around the B(N) atoms, 
and the variation of the exterior angle around the B atoms has a stronger effect on the 
increase in the polarization than that around the N atoms. Thus, the polarization becomes 
totally stronger as the exterior angle increases. 

We next explore the polarization properties with respect to the BN bond length. The 
electronic polarization of the graphitic sheet as a function of the ratio of the cross section 
of the sheets, S/S , is shown in Fig. [51 Here, S = tL x ( y ), t is the effective thickness of the 
sheet j^l, and L x ^ is the length of the supercell in the x(y) direction. In addition, Sq is 
the cross section of the sheet in which the lengths of all the BN bonds are a^. The length 
of the supercell in the direction perpendicular to the external electric field is varied, while 
that in the direction parallel to the field is fixed and the graphitic sheet is kept flat. For 
both directions of the external electric field, the polarization proportionally becomes large 
with increasing lateral length L x ^ . 

To explore the relation between the polarizability and the geometrical deformation due 
to the wrapping into the tubular form, we plot in Fig. [6] the variations of the exterior angle 
of the BN bonds and the cross section of the BN NTs as a function of the diameter. When 
the BN graphitic sheet is wrapped into the tubular form, the exterior angle increases and 
the cross section becomes small with the decrease in the diameter. The external angle of the 
zigzag NTs varies more than that of the armchair NTs; the large variation of the exterior 
angle increases the polarization of the zigzag NTs. In contrast, the decrease in the cross 
section results in the weak polarization of thin armchair NTs. 

IV. CONCLUSION 

We have studied the electronic polarization of BN and A1N NTs and their graphitic 
sheets under finite external electric fields from first principles. We found that the electronic 
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polarization per atom of the zigzag NTs increases as their diameter decreases. On the 
contrary, the polarization of the armchair NTs decreases as the NTs become thin. The 
variation of the polarization with the diameter is associated with specific features of the bond 
arrangement, such as the bond angles around the B or Al atoms and the bond length. Our 
results indicate the possibility of tuning the dielectric property of III-V NTs by controlling 
their diameter and chirality. 
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TABLE I: Polarization per atom (a.u./atom) of graphitic sheets. 





BN A1N 


Polarization per atom 


0.0362 0.0523 
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FIG. 1: Computational models of (a) zigzag NT, (b) armchair NT, and (c) graphitic sheet. Light 
and dark spheres indicate B and N atoms, respectively. E is an external electric field. 

TABLE II: Radial buddings (A) of III-V nanotubes. 





BN 


A1N 


(6,0) 


0.090 


0.119 


(9,0) 


0.057 


0.083 


(12,0) 


0.042 


0.062 


(4,4) 


0.080 


0.113 


(5,5) 


0.061 


0.091 


(6,6) 


0.050 


0.074 
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FIG. 2: Electronic polarization per atom as a function of diameter. The zero of the polarization is 
set to be those of the graphitic sheets, which correspond to 36.2 xlO -3 a. u. /atom for BN NTs and 
52.3xl0~ 3 a.u./atom for A1N NTs. 
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FIG. 3: Computational models of graphitic sheet in which B or N atoms are displaced along the 
direction perpendicular to the graphitic sheet, (i) B atoms are shifted, (ii) N atoms are shifted. 
Gray and dark-blue spheres indicate B and N atoms, respectively. In the top views, B(N) atoms 
are denoted by larger and smaller spheres according to their distance from N(B) atoms. 
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FIG. 4: Variation of electronic polarization, AP = P(#) — P(0), of BN graphitic sheet as a function 
of exterior angle 9 in Fig. [3l (a) and (b) The external electric field is applied along the x direction, 
(c) and (d) The external electric field is applied along the y direction. The curves (a) and (c)[(b) 
and (d)] show the polarizations of model (i) [(h)] in Fig. [3l 




FIG. 5: Variation of electronic polarization, AP = P(S) — P(Sq), of BN graphitic sheet as a 
function of variation of cross section. The computational model is the same as that in Fig. [T] (c). 
The curve (a) [(b)] shows the polarization of the case that the external electric field is applied along 
the x[y] direction and the length of the supercell in the y[x] direction is varied. 
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FIG. 6: Variations of exterior angle of the BN bonds and cross section of the BN NTs as a function 
of diameter of NTs. The solid and dotted curves are the variations of the zigzag and armchair 
NTs, respectively. 
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